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ABSTRACT 

Several properties of the standard a-disc model for active galactic nuclei (AGN) are not 
entirely consistent with AGN observations. As well as such discrepencies, observations 
show evidence for the existence of high mass outflow winds originating from the vicinity 
of the active black hole (BH) . Such winds may originate from various parts of the disc 
and could change the local accretion rate which should alter the emitted spectral 
energy distribution (SED) and affect the global disc luminosity and the BH growth 
rate. The new calculations presented here show the effects of several types of winds 
on the observed and inferred disc properties. Some wind profiles can have a profound 
effect on the observed SED and can perhaps explain the poorly understood deviations 
of AGN spectra from standard disc spectra. We show a factor ~ 2 possible error 
in estimating the disc luminosity and larger deviations in estimating L/L c dd- The 
BH growth rate computed without taking the effects of wind into account may be 
significantly over-estimated. We also suggest a practical way to use the observed SED 
in order to make first order corrections to BH growth rate and account for the effects 
of disc winds. 

Key words: (galaxies:) quasars: general; galaxies: nuclei; galaxies: active; accretion, 
accretion discs 



1 INTRODUCTION 

Several properties of the spectral energy distribution (SED) 
of Active Galactic Nulcei (AGN) can be explained by mass 
accretion via a disc onto a central black hole (BH). AGN 
accretion discs (AD) may vary in geometry, optical depth, 
accretion rate and other parameters. The basic theory of op- 
tically thick, geometrically thin ADs (k nown as a-discs) is 
described in IShakura fc Sunvae v (1973). More recent stud- 
ies that include Comptonization in the disc atmosphere and 
improved trea t ment of vertical structure are described in 
iHubenv et alj (|200lf ) and references therein. Observational 
evidence, mainly in the optical-UV, suggests that the a- 
disc model spectrum is a good approximation of real AGN 
continua. However, there are a number of characteristics 
which cannot be explained by this model. In particular, for 
AGNs with Mbh « lO 8 A/ and M/M edd « 0.3, where 
M is the accretion rate and M e dd is the Eddington ac- 
cretion rate, the model predicts a power-law, F v oc v a , 
with a = 1/3 throughout much of th e optical-NUV con - 
tinuum. IVanden Berk et all (|200lh and iDavis et all (|2007m 
found that for many objects observed in the Sloan Digital 
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Sky Survey (SDSS), the 1350 - 4200A power-law is consid- 
erably softer, with a « —0.5. Observations of t he FUV and 
EUV continua show much softer power- laws (|Telfer et al.l 
120021 ; IScott et al.l 120041 : IShull et al.|[2012h . Apart from the 
Q-disc continuum, the overall SED is comprised of contribu- 
tions from a number of other components incl uding emission 
by ho t dust in the NIR part of the spectrum (|Mor fc Netzej 
l201lf ). a conglomeration of hundreds of broad Fell lines 
and Ba l mer continuum emission (the "Small Blue Bump". 
iGrandil Il982l; iNetzer et al.l Il985l) and a power-law x-ray 
source (|Laor et al.lll997l ). |jin et al.l (|201ll ) tried to combine 
these processes in order to model a large sample of SDSS 
spectra but found that for a great number of objects the un- 
derlying continuum was softer than they were able to m odel. 
Similar conclusions were reached bv lDavis et alj (j2007) , who 
attempted to compare observed spectral properties of a large 
sample of AGNs to those predicted by their models. A de- 
tailed overview of the variat ions between AGN theor y and 
observations can be found in iKoratkar fc Blaesl (|l999h . 

AGN ADs are also probably associated with high ve- 
locity winds. X-ray observational evidence f or such winds 
is giv en in numerous public at ions including Pounds et al.l 
d2003h. iReeves et al l (|2003l h iTombesi et al.l and 
ITombesi et all (|201 ll) and references therein. These stud- 
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ies claim to detect winds, based on Fe K-shell transitions 
with velocities reaching considerable fractions of the speed 
of light. Broad absorption line (BAL) outflows provide more 
evidence for materia l which may be ejected from the vicin- 
ity of the BH (e.g. ICapellupo et al l l201ll . and references 
therein). Statistical analysis of AGNs shows a correspon- 
dence between object s exibiting BAL outfl ows and softened 
optical-UV spectra (| Ganguly et al.l 120071 ). Slower winds, 
which probably originate further from t he BH, may be de- 
tected through ion or m olecular lines fe.g. lStorchi-BergmannI 
l2010l ; ISturm et al.l201l] , and references therein) . iKing! (|2010l ) 
showed that very general considerations predict high ve- 
locity winds expected from AGNs. All these findings sug- 
gest that winds are present in a considerable fraction of 
AGNs. The origin of these winds, especially those with 
high velocities, could be from regions close to the BH 
and are po ssibly connected with, or ejected from, the AD 
itself fe.g IKing fc Pounds! 120031 ; iProga fc Kallmanl 12004 
ISim et al.ll2008l ). 

Winds from AGN ADs may be caused by a num- 
ber of physical processes. Three main scenarios have been 
suggested: radiation driven, thermally driven or magneti- 
cally driven winds. Radiat ion driven winds have been stud- 
ied in numerous pap ers l|Shieldsl Il977l : IProga et all [2000; 
IRisaliti fc ElvisI |2010| . and references therein). Such pro- 
cesses accelerate gas through radiation force on resonant 
line transitions and are able to produce outflows with ve- 
locities of the order ~ 10 4 km ■ s -1 a rising fro m radii within 
the in ner AD. lBegelman et al.l (|l983l ) and later I Woods et all 
( 1996) proposed the possibility of thermal winds caused by 
heating of the disc by a central x-ray source driving ther- 
mal outflows at large radii. Many other studies of such 
wind s also exist jKrolik fc Vrti lek 1984; Che louche fc Netzerl 
l2005l : lEverett fc Murravll2007i . and references therein). Mag- 
netically driven winds seem to be the best candidate for 
the origin of high velocity winds emitted from the in- 
ner radii of ADs. The details o f such systems have been 



described by many researcher s (Blandford & Pavnel 



Contopoulos & Lovelace 1994; Konigl & Kartic 1994 



1982; 



and 

references therein). A common result is a self similar wind 
with a constant mass los s over each decade of the disc ra- 
dius. iBlandford fc Payne! (|l982l ) showed that for standard 
AGNs, only a small fraction of the disc mass is ejected by 
such winds. 

The standard a-disc model assumes a constant accre- 
tion rate, M, throughout the disc. Significant outflows must 
alter the local accretion rate, in turn altering the energy flux 
from different radii. This will change the emitted SED and 
will also affect the simple relationship between disc lumi- 
nosity, L, and the BH growth rate, Mbh- In this paper we 
investigate two major aspects of the possible effects of disc 
outflows on AGNs: the direct effects on the emitted SED, 
and the ability to estimate BH growth rate based on mea- 
surements of L and the global SED. The structure of the 
paper is as follows: In §2 we give a detailed description of 
our numerical model which includes the effects of disc winds. 
In §3 we give the results of our numerical models and show a 
number of examples which demonstrate trends in the SED, 
in L and in the BH growth rate as functions of various AGN 
parameters and wind profiles. Finally, in §4 we discuss our 
results and compare them to AGN observations. 



2 CALCULATIONS 



2.1 Standard Q-Disc Models 



Our calculations are similar to several earlier AD calcula- 
tions based on the standard a-disc model, with the addition 
of GR corrections and treatment of Comptonization in the 
disc atmosphere dBardeen et al. 19721: Laor fc Netzerl 1 19891 : 



IHubenv et al1l200ll : IDavis fc Laorll201ll ~ For completeness 
we summarise the ingredients of the model that are needed 
for the new calculations. 

Three parameters define the emitted SED for a stan- 
dard a-disc: BH mass (Mbh), accretion rate (M) and nor- 
malized BH spin (a). Viscous torque m ay be parameterized 
via th e a V i ac parameter as defined in IShakura fc SunvaevI 
(| 19731 ). The standard model assumes a constant M, inde- 
pendent of radius (R). The BH spin defines the innermost 
stable circular orbit (ISCO) around the central BH and sets 
the efficiency, r/, of the conversion of rest energy to electro- 
magnetic radiation (L = r/Mc 2 ). Angular momentum within 
the disc is transferred outward causing mass accretion onto 
the central BH and viscous dissipation leads to heating of 
the gas and emission of radiation. A local balance is attained 
between work done by torque forces between adjacent differ- 
ential annuli, gravitational potential energy, kinetic energy 
and emitted radiation. The torque forces, which are powered 
by gravitational energy, may correspond to either positive or 
negative work done by the torque. For the standard a-disc, 
the work done by torque changes sign at a radius R = 9/4 
R g , where R g = GMbh/c 2 is the gravitational radius. Above 
this radius, energy is contributed to the emitted radiation 
and below this radius, energy is removed. 

A derivation of the emissivity associated with viscous 
dissipation may be done using the Navi er-Stokes equations 
for a non-compressable viscous fluid (see lClarke fc Carswelll 
l2003h. We derive the same result as theirs using the approach 
of lFrank et~aT] (|2002l ). The emmisivity of the disc at various 
radii depends on the local gravitational energy and work 
done by the local torque. We show below that both these 
components are altered when a wind is taken into account. 
The total emmisivity is 



where 



D(R) = D N (R) + D G (R), 
D G {R) = M- GMbh 



(1) 



(2) 



871-i? 3 

is the emmisivity component due to loss of gravitational 
and kinetic energy. Dn(R) is obtained by considering the 
luminosity due to the work done by the torque, 

Ln{R) = -(^E] 1/2 f R m (GMbh) 1 ' 2 d (r 1 ^ 2 



\ R 3 J Ju in 

(3) 

where the integrand is the differential of the local torque. 
The emmisivity associated with this component is then 



D N (R) = 



(4) 



AirRdR 

The effective temperature of the disc as a function of radius 
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■ M 8 = 1 M/Mdd = 0.3 a = 0.998 

■ M 8 = 0.1 M/Metd = 0.3 a = 0.998 




Figure 1. Calculated SEDs for a-discs with various values of Mg, M /M edd an d a, and comparison to earlier models. Blue and red points 
correspond to values from Hubeny et al. (2001) for the same parameters as used in our models (blue and red lines respectively). They 
fit well with our models for frequencies below ~ 3 X 10 16 Hz and are harder than our models for higher frequencies (see text). 



-3/4 



(5) 



where r = R/R g and f(r) is a function of order unity (for 
large values of r) which depends on the GR terms. Assuming 
local blackbody (BB) emission, the isotropic L v is 



L v = / 4irR(nB v (T))dR= 



RdR 



R exp(x) — 1 ' 
(6) 

where x — hu/kT and Ri„ and R ou t are the innermost and 
outermost radii of the disc respectively. The dependence of 
L v on temperature and disc geometry is isolated in the in- 
tegral 

^ RdR 

jR in exp(x) - 1 

Our AD calculations include Comptonization of the 
emitted radiation in the AD atmosp here. We model Comp- 
tonization in the manner applied bv lLaor fc Netzerl |l989), 
namely, we calculate the parameter y (their eqn. 29) which 
is a function of r and v and is dependant on the viscosity 
parameter, a v i 3C , and on the gas c omposition. W e have com- 
pared our values of y to those of iHubenv et al.l (120011 ) and 
find them to be in good agreement. 



Fig. [T] shows some examples of our standard disc 
calculations and a comparison with those calculated by 
IHubenv et alj (|200lh . The agreement is good, except for 
very high frequencies where our SEDs are somewhat softer, 
mainly because we do not account of radiative transfer 
within the disc atmosphere. Fig. [T] also demonstrates the 
known dependence of the SED shape on Mg (defined by 
M s = M B h/(1O 8 M )), M/M edd and a. BH spin governs the 
radius of the ISCO, changing the spectral shape through the 
boundaries of the integral 6{u) (eqn.[7|). Mbh and M change 
the spectral shape through T(r). 

The a-disc model leads to a simple relation between the 
luminosity and M, 

-Lodd M edd 

where L — L v dv, L c dd scales linearly with Mbh and 
M^d depends on Mbh and a. As shown below, this simple 
relation breaks down when M is no longer constant. 



2.2 Wind Ejecting a-Disc Models 

We now explore how winds ejected directly from the sur- 
face of the disc can affect the SED. Ejection of mass from 
various radii causes M to change as a function of r, thus 
M is replaced by M (r). We define M edd exactly as above, 
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• wind type 1 (self similar') 
■ wind type 2 (step function) 

wind type 3 (linear) 

M=const 



300 



400 



500 



R/Rg 



Figure 2. Example M (r) profiles for models with Mg = 1, M in /M edd = 0.1, M out /M edd 
50% of total wind mass has been ejected from the disc. All models shown here have T\/2 '■ 
line is a standard a-disc with constant M/M edd = 0.1. 



0.3 and a = 0.9. is the radius at which 
50 (vertical dashed line). The dotted black 



thus, it's value is independent of the existence of wind. We 
assume the local M (r) to be approximately constant within 
each differential annulus, allowing us to calculate the local 
emmisivity of each annulus in a way similar to the stan- 
dard a-disc model. We assume a certain mass inflow rate, 
Min, at the ISCO. We then specify a radius dependent mass 
outflow rate (hereafter "wind") which defines M (r). We ex- 
perimented with three functional forms for the wind profile: 
(1) a self similar w i nd fu nction, such as that described by 
iBlandford fc Pavnd <|l982r ). (2) a step function wind profile 
and (3) a wind function which decreases linearly with radius. 

We define ri/2 to be the radius at which 50% of total 
wind mass has been ejected from the disc (see fig. The 
self similar, type (1) wind function, is of the form 



dM wind (r) 
dr 



r 



(9) 



where ai is a normalization parameter and depends on Mi n 
and M ov .t- For such a wind and for Ms = 1, a = 0.9 and 
a disc with a maximum radius of 1080 R g , we get r 1( / 2 = 
50. The type (2) wind profile resembles a step function and 
depends on two parameters, fei and i>2, 

dM "7 d(r) = h (1 - [1 + exp (-26 2 (r - r m ))r') , (10) 
dr 

where r m is a parameter which defines rx/i- The type (3) 



function is a linear function of the form 

dM wind (r) 



dr 



(ci - c 2 • r), 



(11) 



where ci and C2 are chosen such that r^/i has the same 
definition as above. In all three cases, the wind function 
defines M (r) according to 



M (r) = M in + 



dM u 



i{r). 



(12) 



Diagrams of M (r) for some of the examples discussed in this 
paper are shown in fig. [2] We found very similar results for 
winds of types (2) and (3) (see details in $3]), and significant 
differences for winds of type (1) . In this work we show results 
for wind types (1) and (2). 

The varying M (r) has a number of effects on the disc 
geometry and the emitted SEEQ. The removal of energy from 
the AD changes the balance between the gravitational po- 
tential and torque energies within the disc. Dc(r), Djv(r) 
and T(r) may be calculated locally according to equations 
[2][3]|4]and[5l with the values of M (r) at each annulus mod- 
ified according to the wind function. The gravitational com- 
ponent, Dg{t), is simply proportional to M (r). The depen- 
dence of the torque component on M (r) is less straightfor- 

Note that we do not consider the fate of the wind after it leaves 
the disc. 
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Figure 3. Plots of emmisivity for various type (2) wind functions for models without Comptonization and GR effects. All models have 
Mg = 1, Mi n / M e dd = 0.1, a = 0, a v isc = 0.1, f>2 = 0.4 and = 10. Blue lines correspond to a standard a-disc model with constant M. 
Solid lines correspond to emissivity associated with gravitational and kinetic energy loss, Dq(t). Dashed lines correspond to emissivity 
associated with work done by torque, Dpf(r). Dq(t) is modified locally proportionally to M (r) at each radius. Regarding Djy(r), the 
wind profile changes M(r) such that the integral in eqn. \3\ is modified, thereby changing the form of D^(r). As more mass is lost 
throughout the disc, torque forces do negative work and consume energy over a larger range of radii (between r'i n and the radius where 
D N (r) =0). 



ward. The radius where Dn(t) = is no longer r — 9/4. Fig. 
[3] shows Djv(r) and Dc{r) for four different values of M ou t 
and the same M in , calculated without including the effects 
of Comptonization and GR (discussed below). The modifi- 
cations of the emissivity change T(r) , and accordingly L„ is 
altered through the modification of the function 6(u) (eqn. 
0. 

To the best of our knowledge, GR corrections, which 
have the effect of cooling the inner radii of the disc, have 
never been calculated for a wind ejecting a-disc. We have 
therefore adopted the following procedure. We calculated 
the function 



using the standard a-disc model (i.e. M = const), where 
TNewt(r) and T re i(r) are Newtonian and GR temperatures 
respectively. The values of the function F(r) depend only 
on the gravitional potential and therefore only on a. For the 
wind calculations we obtain T(r) numerically and multiply 
by F(r). This has the effect of cooling the inner radii of our 
models, where relitavistic effects are most prominent. 

In wind ejecting discs, the BH growth rate no longer 



depends on the accretion rate at the outer radii of the disc, 
only on Mj„. Thus, eqn. [8] is no longer valid and measure- 
ment of L is no longer a direct measure of Mbh even if a is 
known. We later discuss the implications and possible errors 
in estimating BH growth rates when the effects of winds are 
ignored. 

We tested our model for a large range of A/bh, a, M (r) 
and a V i SC . M (r) is characterized by Mj n , M out , ri/a and 
the wind function. We show here only the most interesting 
results and only those with wind function types (1) and 
(2) (eqns. [9] and 1 10 p . The results are almost independent 
of our choice of a V i SC which we have chosen to be 0.1 in all 
our calculations. In table Q] we give details of the models 
shown in this paper. The models have all been calculated 
with Mg = 1 and each with differing BH spin, Mi n , Mout and 
values of ri/2- It is convenient to compare different models 
by deviding them into three catagories: (A) keeping M;„ and 
Mout constant while differing r lt / 2 , (B) keeping Mi„ and ri/2 
constant while differing M m t an d (C) keeping M out and rw2 
constant while differing Mi„. Models in table[T]are classified 
according to BH spin, category (A) or (B), and wind type 
(1) or (2). Comparing models 1.3 and 1.8, models 2.3 and 
2.8 and models 3.3 and 3.8 gives category (C). For models 
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Table 1. Parameters of the models shown in this paper. All models have Ms = 1 and a v i ac = 0.1. Models with wind type (2) fegn. [T0) l 
have 62 = 0.4. Models 1.1, 2.1, 1.5 and 2.5 are standard a-discs with constant M. All other models are wind ejecting a-discs. Models 1.2 
- 1.8 and 2.2 - 2.8 were calculated with wind type (2), models 3.3, 3.6, 3.7 and 3.8 were calculated with wind type (1) and have the same 
parameters as models 2.3, 2.6, 2.7 and 2.8 respectively. Models with wind type (2) are divided into categories with either the same Mi n 
and M ut an d differing rj/ 2 (category A) or the same M; n and rj/ 2 and differing M out (category B). Category (C) models are 1.3, 1.8, 
2.3, 2.8, 3.3 and 3.8 which have the same M ou t and Tin and differing Mi n . 



Model 


BH spin (a) 


M in 


M out 


r l/2 


Wind Type 


Category 






[M edd ] 


[M edd ] 






Model 1.1 





0.1 


0.1 








Model 1.2 





0.1 


0.3 


10 


2 


A 


Model 1.3 





0.1 


0.3 


50 


2 




Model 1.4 





0.1 


0.3 


250 


2 




Model 1.5 





0.01 


0.01 








Model 1.6 





0.01 


0.03 


50 


2 


B 


Model 1.7 





0.01 


0.1 


50 


2 




Model 1.8 





0.01 


0.3 


50 


2 




Model 2.1 


0.9 


0.1 


0.1 








Model 2.2 


0.9 


0.1 


0.3 


10 


2 


A 


Model 2.3 


0.9 


0.1 


0.3 


50 


2 




Model 2.4 


0.9 


0.1 


0.3 


250 


2 




Model 2.5 


0.9 


0.01 


0.01 








Model 2.6 


0.9 


0.01 


0.03 


50 


2 


B 


Model 2.7 


0.9 


0.01 


0.1 


50 


2 




Model 2.8 


0.9 


0.01 


0.3 


50 


2 




Model 3.3 


0.9 


0.1 


0.3 


50 


1 




Model 3.6 


0.9 


0.01 


0.03 


50 


1 




Model 3.7 


0.9 


0.01 


0.1 


50 


1 




Model 3.8 


0.9 


0.01 


0.3 


50 


1 





with wind type (1), the radius of the outer edge of the disc 
was chosen such that ri/2 = 50. Note that models with 
the same Mi„/ M edd and different a (e.g. models 1.1 and 
2.1) do not have the same accretion rate, M, because M edd 
depends on spin. Each model was used to calculate a number 
of parameters which define the shape of the SED, the total 
luminosity, L, and the BH growth rate, Mbh = (1 — fj)Mm- 

Energy conservation requires that the total loss of grav- 
itational and kinetic energy within the disc is equal to the 
total emitted luminosity. Our calculations do not consider 
the fate of the mass within the wind after it has left the disk 
and hence it only includes half the energy required for the 
wind mass to escape the BH potential. There are a num- 
ber of possibilities regarding the fate of this mass which are 
consistent with our model. One such possibility is that the 
ejected wind rises above the surface of the AD, but does not 
leave the physical vicinity of the BH and therefore does not 
require additional kinetic energy. This wind may re-enter the 
disc at larger radii or may accumulate and become ionized in 
regions above the disc (creating a "failed wind"). Another 
possibility may be that the wind is not isotropic (as sug- 
gested by many models) and that it escapes the system by 
absorbing radiation propogating at angles closely parallel to 
the disc. Therefore, if we view the AD at angles close to face 
on, the radiation we observe may not be influenced by this 
effect. Of course, there are other possibilities that involve 
additional acceleration of the wind by part of the radiation 
that we included in our luminosity budget. A detailed dis- 



cussion of these and similar possibilities is beyond the scope 
of our paper. 



3 RESULTS 

In the following section we analyze each model's SED and 
compare the spectra of discs with various parameters to one 
another. Two major aspects of the wind ejecting disc models 
are discussed: the spectral changes, and the effect of mass 
outflow on the ability to evaluate Mbh- Results of the spec- 
tral changes are expressed by showing the approximated 
power-law behaviour, F v tx v a , within various wavelength 
ranges and by the mean weighted ionizing photon energy of 
the Lyman continuum, < Eionphot >• The degree of inaccu- 
racy in evaluation of Mbh from the observations depends 
on our understanding of the effect of winds on the SED and 
on the ratio of the measured L to the luminosity that would 
have been observed from a standard (constant M) a-disc 
with M = Mi n . We denote this luminosity Lo, which for a 
wind ejecting model satisfies the equation 

J±. = J^IL . (14) 

Lodd Medd 
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Table 2. Continuum power-laws (L v oc v a ) and mean ionizing photon energies for models shown in table [T] 



Model 


a 


a 


a 


a 


^ ^ionphot -- > 




456— 912A 


912 — 1450/1 


1450 — 4200/1 


4200 — 5100A 














[eV] 


Model 1.1 


-1.52 


-0.53 


-0.044 


0.19 


21.32 


Model 1.2 


-1.31 


-0.43 


0.016 


0.33 


22.07 


Model 1.3 


-1.66 


-0.75 


-0.21 


0.20 


21.10 


Model 1.4 


-1.55 


-0.60 


-0.20 


0.063 


21.27 


Model 1.5 


-3.58 


-1.47 


-0.40 


-0.0053 


17.86 


Model 1.6 


-3.62 


-1.61 


-0.61 


-0.17 


17.83 


Model 1.7 


-3.65 


-1.80 


-0.71 


-0.15 


17.81 


Model 1.8 


-3.41 


-1.66 


-0.56 


0.060 


18.06 


Model 2.1 


-0.61 


-0.14 


0.097 


0.21 


26.47 


Model 2.2 


-0.79 


-0.26 


0.038 


0.22 


25.53 


Model 2.3 


-0.72 


-0.35 


-0.13 


0.085 


23.098 


Model 2.4 


-0.63 


-0.19 


-0.031 


0.025 


26.40 


Model 2.5 


-1.59 


-0.58 


-0.086 


0.11 


21.09 


Model 2.6 


-1.64 


-0.70 


-0.31 


-0.12 


21.01 


Model 2.7 


-1.80 


-1.00 


-0.57 


-0.23 


20.78 


Model 2.8 


-2.03 


-1.26 


-0.61 


-0.13 


20.44 


Model 3.3 


-0.69 


-0.24 


-0.012 


0.12 


26.14 


Model 3.6 


-1.63 


-0.67 


-0.19 


-0.0030 


21.068 


Model 3.7 


-1.61 


-0.72 


-0.26 


-0.064 


21.24 


Model 3.8 


-1.41 


-0.65 


-0.24 


-0.022 


22.041 



3.1 Effects of Wind on the Spectral Energy 
Distribution 

We examine the SED shape by fitting a power-law, F v oc 
u a , to four wavelength bands: 456 - 912A, 912 - 1450A, 
1450-4200A and 4200-5100A. For each model we calculate 
these four values of a and the mean ionizing photon energy 
for photons with energies above 13.6eV. SEDs for several 
models are plotted in figs. [4] and [5] The results for a and 
the mean ionizing photon energies are listed in table [2] 

For a wind ejecting disc with given BH mass and spin, 
the spectral shape is influenced by M (r) which can be char- 
acterized by the values of Mm, M out and r ± / 2 - In general, a 
large accretion rate throughout much of the disc has the ef- 
fect of causing higher total luminosity and shifting the peak 
luminosity toward higher frequencies. The closer the high 
accretion rate reaches to n n , i.e. the lower the value of ri/2, 
the more high energy UV radiation will be emitted. The low 
energy, optical-NIR radiation, is emitted mainly from large 
radii and depends less on riy 2 . 

The SEDs of the category (A) models, 1.1 - 1.4 and 
2.1 - 2.4, are plotted in fig. [I] All have M in /M edd = 0.1 
and M out / M e dd = 0.3 (except for standard Q-discs 1.1 and 
2.1 which have constant accretion rates). The shapes of the 
SEDs of wind ejecting models as compared to standard mod- 
els 1.1 and 2.1 are influenced by two competing effects: 
higher total luminosity together with a shift of the peak 
luminosity bluewards, and a sharper fall in L v at high fre- 
quencies, caused by reduction in accretion rate at small radii. 
These two effects tend to cancel each other out in their in- 
fluence on the shape of the UV spectrum as can be seen by 
the similar values of mean ionizing photon energies in the 
last column of table[2] Comparing spectral shapes of models 
with wind to those of models without wind we see that the 



same two effects create very similar spectral shapes between 
models 1.1 and 1.4 and between models 2.1 and 2.4. As r 1( / 2 
grows, the spectral shape changes most in the NIR-optical 
continuum with the UV continuum less influenced because 
of less accretion at small radii. 

Models compared from category (B) are 1.5 - 1.8 and 
2.5 - 2.8, all with M in /M e dd = 0.01 and r 1/2 = 50 (except 
for standard a-discs 1.5 and 2.5 for which ri/2 has no mean- 
ing). Their SEDs are plotted in fig.[S] The growing M ou t has 
it's greatest effect on amplification of the NIR-optical con- 
tinuum and a lesser effect on the UV continuum. Spectral 
changes between models with wind to those without wind 
are more prominent for BHs with a = 0.9 than for BHs with 
a — (as can be seen, for example, by comparing model 
1.5 to 1.8 and model 2.5 to 2.8). This is because for higher 
spin values, smaller radii are influenced by the wind. The 
changes are most prominent at lower frequencies. 

Comparison of models from category (C) (models 1.3, 
1.8, 2.3, 2.8, 3.3 and 3.8) have different behaviour. For such 
AGNs, with the same M out and r\/2 and different values of 
Mi n , the most prominent spectral changes are in the UV 
spectrum. The reason for this is the very different accretion 
rates at small radii and the very similar accretion rates at 
large radii. We will later show that this characteristic of 
wind ejecting discs enables us to best constrain BH growth 
rate for AGNs with known luminosities. 

We find that our wind models are able to soften the 
optical-UV SED (decrease the value of a) creating a result 
more similar to observations of l arge samples of AGN s . This 
can be seen by comparing the IVanden Berk et al.l (|200ll ) 
composite to our models in fig. [5] Changes in the observed 
a can also be caused by a different combination of BH mass 
and accretion rates, in cases without wind. We have there- 
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Figure 4. Comparison of SEDs for models of category (A). All have the same M; n = 0.1 and therefore the same Mbh- Left: SED 
comparison for models 1.1 - 1.4 with a = 0, Right: SED comparison for models 2.1 - 2.4 with a = 0.9. Solid blue lines - models 1.1 and 
2.1; dashed red lines - models 1.2 and 2.2; dotted dashed green lines - models 1.3 and 2.3; dotted dashed magenta lines - models 1.4 and 
2.4. Vertical dotted lines mark boundary wavelengths used for calculations of power-laws in table [2] 




Figure 5. Comparison of SEDs for models of category (B). All have the same Mi„ = 0.01 and therefore the same Mbh- Left: SED 
comparison for models 1.5 - 1.8 with a = 0, Right: SED comparison for models 2.5 - 2.8 with a = 0.9. Solid blue lines - models 1.5 and 
2.5; dashed red lines - models 1.6 and 2.6; dotted dashed green lines - models 1.7 and 2.7; dotted dashed magenta lines - models 1.8 and 
2.8. Vertical dotted lines mark boundary wavelengths used for calculations of power-laws in table [2] Also shown in black is a composite 
SED from Vanden Berk (2001) with arbitrary normalization. Note that the reddening of the IR-optical continuum of the composite is 
caused by a significant contribution of the host galaxy. 



fore searched for standard a-disc models which give similar 
SEDs to those of our modified model. We allowed freedom 
of all relevant disc parameters (Mbh, M, a) and created 
SEDs with optical-UV power-laws a ° ~ —0.5. Fie. 

1450-4200A ° 

[6] shows contours of a » = —0.5 for a large range of 

1 1450-4200A ° ° 



Ms, M and various values of a. We were only able to model 
such soft optical and UV spectra by choosing Ms > 1 or 
M / M e dd -C 0.1 which both have the effect of softening the 
optical continuum. We were unable to model soft enough 
optical-UV spectra with M$ « 1 and M/M edd w 0.1 - 0.3. 
This has important implications, since such soft continua 
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Figure 6. Contours of a 145() 42 ooA = ~~ "'^ ^ or ^a^dard c*-discs with varying Ms and M /M e dd ar e shown here. Each contour corre- 
sponds to a different value of BH spin, a. The grey area corresponds to the parameter space where the standard a-disc cannot model 

a • = —0.5 for a ^ 0. Low accretion rates and large BH masses soften the optical-UV continuum. Standard a-discs with M& 1 

1450-4200A 6 1 s 

and M/M e dd ~ 0.3 cannot create c* 1450 42 ooA = — ^ so ^ °P^ ca l power-law corresponding to observations can be achieved by the 
standard a-disc model only for large Mbh ° r l° w accretion rates. 



are observed in objects where the BH mass is observed to 
be around Mg = 1. 

3.2 Effects of Wind on the Estimation of Black 
Hole Growth Rate 

As discussed above, for a wind ejecting disc, there is no sim- 
ple way of calculating Mbh by measuring L. Since Mbh 
depends only on spin and Mi n , and since models in each 
category have the same Mi n , all those with equal spin have 
the same BH growth rate. However the total L and the spec- 
tral shape of each can greatly vary. Table[3]shows a summary 
of factors influencing the correlation between Mbh and L. 
L/Lq is the ratio of the luminosity of the model to that of a 
standard a-disc model with M = Mi n (Lo is defined in eqn. 
I14[) . For models with the same Mi n and M out but different 
r-i/2, a disc which loses most of its mass at smaller radii 
(models 1.2 and 2.2) result in larger L/Lo than a disc which 
loses most mass at larger radii (model 1.4 and 2.4), because 
the luminosity of the latter case is smaller. Obviously, L/Lo 
can reach extremely large values when M ou t/Mi„ is large. 
For our most extreme models, L/Lo ~ 7. 

The standard way to estimate L from observations is 
to use a wavelength dependent bolometric correction fac- 



tor, b(X) — Lj\L\ at a chosen wavelength. For example 
6(5100A) is used to estimate L in low redshift AGNs. Our 
calculated values of b(5100A) and the standard b olometric 
correction calculated by the M arconi et alj |2004l ) method 
for 5100A (6m) are given in table [3] We note that the 
iMarconi et al.l |2004l ') method includes contributions to L 
from both the optical-UV and x-ray continua, whereas we 
do not include the x-ray continuum in our models and there- 
fore calculate L without the x-ray contribution. The same 
method, when used for L(x — ray) , integrated over the x-ray 
continuum, results in Lix — ray)/L in the range ~ 1 — 30% 
depending on L and on the energy range used for integration. 
Therefore, the fact that we do not take this luminosity con- 
tribution into account may lead to differences in bolometric 
correction factor of up to ~ 20% for extreme cases and has 
a minor effect on the bolometric correction calculations for 
most of our models. The ratio 6a//6(5100A) depends on the 
spectral shape of the AGN. Values of this ratio close to unity 
imply accurate estimates of the disc luminosity even in the 
presence of wind. Multiplying L/Lo by &m/6(5100 j 4) (3rd 
and 6th columns of the table) gives the ratio (6mALa)/£o 
which is the error in evaluation of Mbh if only a measure- 
ment of AL\ is available. These values are given in the last 
column of the table. 
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Table 3. Changes in factors effecting estimation of Mbh f° r wind emitting a-disc models. Lq is defined in eqn. If 41 fe(5100A) is the 
bolometric correction calculated according to b(X) = L/\L\\ bM is the bolometric correction as calculated by the Marconi (2004) method 
at 5f00A; 6ji//6(5100A) is equal to (bju \L\)/ L; bM \L\/Lq is the error in evaluation of Mbh if only a measurement of \L\ is available 
and therefore may be written M B H eatimated / M B Hreal- 



Model 




L/L 


6(5100A) 




b M /b(5100A) 


but^Lx/Lo 


Model f.f 


0.1 


1 


8.95 


7.41 


0.83 


0.83 


Model 1.2 


0.23 


2.37 


10.41 


7.028 


0.68 


1.60 


Mftfjpl 1 Q 

iviooei i.o 


14 


1 42 


7 94 


7 10 


u.yo 


1 40 


Model 1.4 


0.11 


1.11 


7.50 


7.25 


0.97 


1.07 


Model 1.5 


0.01 


1 


4.69 


8.61 


1.83 


1.83 


Model 1.6 


0.014 


1.42 


3.94 


8.17 


2.08 


2.95 


Model 1.7 


0.029 


2.90 


3.60 


7.61 


2.11 


6.13 


Model 1.8 


0.071 


7.14 


4.07 


7.16 


1.76 


12.56 


Model 2.1 


0.1 


1 


16.50 


7.80 


0.47 


0.47 


Model 2.2 


0.16 


1.66 


13.79 


7.37 


0.53 


0.87 


Model 2.3 


0.12 


1.20 


11.71 


7.46 


0.64 


0.76 


Model 2.4 


0.103 


1.05 


13.79 


7.65 


0.55 


0.58 


Model 2.5 


0.01 


1 


8.41 


9.17 


1.091 


1.091 


Model 2.6 


0.012 


1.20 


6.57 


8.76 


1.33 


1.60 


Model 2.7 


0.019 


1.89 


4.88 


8.12 


1.66 


3.14 


Model 2.8 


0.038 


3.86 


4.48 


7.57 


1.69 


6.52 


Model 3.3 


0.13 


1.35 


13.71 


7.49 


0.55 


0.73 


Model 3.6 


0.013 


1.35 


7.40 


8.76 


1.18 


1.59 


Model 3.7 


0.025 


2.55 


6.84 


8.15 


1.19 


3.043 


Model 3.8 


0.059 


6.00 


7.43 


7.61 


1.02 


6.15 



We see in table [3] that for models with M ou t/Mi„ = 3, 
L/Lq is in the range 1 - 2.37 for models with a — 0, and in 
the range 1 - 1.66 for models with a = 0.9. The minimum 
values correspond to large r\n and the maximum values 
to small ri/2- The values of 6m/6(5100A) are in the range 
between 0.47 - 2.11 for all models, indicating a possible devi- 
ation of factor ~ 2 if wind is not taken into account. Values 
of bmXLx/Lo for models with M out /Mi n = 3 range between 
1.07 - 2.95 for models with a — 0, and between 0.58 - 1.6 for 
models with a — 0.9. In general, higher spin values give lower 
ratios of L/Lo, &a//&(5100A) and their product, correspond- 
ing to smaller errors in estimation of Mbh from L assuming 
a standard (constant M) a-disc. Models with larger ratios 
of M ou t/Mi n may have even higher errors in estimation of 
Mbh- For our most extreme cases, with M out /Mi n = 30, 
we get b M \L\/ Lq « 12.5. 

For a standard a-disc model, w hen Ms is known, on e 
may estimate M using the method of iDavis fc Laorl i|201lf ). 
utilizing the fact that L„ is insensitive to BH spin at long 
wavelengths. For a wind ejecting model, use of this method 
may cause significant errors in estimation of the accretion 
rate at the ISCO and consequently the BH growth rate. 
We have calculate d the r atio of M calculated according to 
the IDavis fc Laorl l|201lf ) method to the real Mi n for all 
our models and find significant deviations for models with 
Mout/Min > 1. 

We have calculated the values of L/L c dd for models 
with Ms = 1, a = 0, M out /M e dd = 0.3 and all wind 
profiles. This results in a continuous range of values for 
M in /M edd = 0.01 - 0.3 and r 1/2 = 10 - 250. The blue solid 
lines in fig. [7| show Mi„/M e dd as a function of L/L e dd for 



two values of rx/2, 10 and 250. Each coordinate within these 
boundaries corresponds to a possible model with a certain 
Mi n and L. For a given L/L c dd, there is a broad range of 
allowed values of Mi n /M e dd for various wind profiles. The 
diagram shows that there is a value of L/L e dd for which 
there is a maximum range of Mi„/ M e dd, corresponding to 
the luminosity for which we have maximum uncertainty in 
the accretion rate profile (i.e. in the form of M (r)). We have 
calculated the values of a for the four wavelength ranges 
shown in table [2] for all models within the parameter space 
confined by the blue lines of fig. [7| (i.e. those with Ms = 1, 
a — and M out /M e dd = 0.3). The figure shows contours of 
q ■ (left) and a ° (right), a is most sensitive 

456-912A v ' 4200-5100A v ° ' 

to changes in Mi n in the UV part of the spectrum. We show 
this result only for M ou t/M e dd = 0.3, but it holds true for 
all values of M ou t, with the actual shapes of the contours 
varying for different values of M out . In each panel of fig. 
[7] we show in dashed black lines the borders of L/L e dd vs 
Mi n /M e dd for AGNs with M 8 = 1 and M ou t/M ed d = 0.2. If 
measurements of L and of a in the UV are available for an 
AGN, one can constrain the value of Mbh by utilizing the 
sensitivity of a in the UV to Mi„ and taking into consider- 
ation the unknown parameters M ou t and a. 



4 DISCUSSION 

The wind ejecting disc model presented in this work is mo- 
tivated by two observed phenomena: discrepancies between 
observed optical-UV SEDs of AGN to those predicat ed by 
standard a-disc models (as found bv lDavis et alj |2007). and 
evidence for high velocity winds which may be ejected from 
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Figure 7. Sensitivity of UV and optical power-laws to Mi„ and L. Shown in blue solid lines are borders of possible values for Mi n / ^ledd 
as a function of L/Ledd for AGNs with Mg = 1, M ou t/ = 0.3 and a = and all possible wind profiles (i.e. all values of rj/ 2 ). 
Each coordinate within these boundaries corresponds to a possible model with a certain M;„ and L and varying M (r) profiles. Shown 
in dashed black lines are the same borders for AGNs with Mg = 1, Mout / M^d = 0.2 and a = 0. Colored contours are shown of constant 
power-law values for two wavelength ranges (UV and optical) from table[2]- Left: ct * and Right: ct These contours are 

4o0 — y ±^.r± 4^UU — □ -LUU-/i 

shown for the M ou t/ M e dd = 0-3 models only, they differ from those of the M ou t/ 1 M e dd = 0.2 models but retain the same general trends. 

ct ■ (UV) is far more sensitive than a ° (optical) to changes in M; n , therefore, knowledge of ct • best constrains 

456-912A v ' 4200-5100A 1 ' & ln ' ' 6 456-912A 

Mi n for a given L. 



the inner AD. The models calculated here are not intended 
to represent real disc winds. The properties of such winds 
may depend, among other factors, on disc turbulence, the 
rate and direction of mass inflow, disc geometry and mag- 
netic fields, all not very well understood. Our limited scope 
calculations are meant to study the influence of general prop- 
erties of winds and their radius dependent mass loss, on the 
fundamental disc properties such as total luminosity, SED 
shape and BH growth rate. Obviously, the calculations do 
not cover all such possibilities. Furthermore, we stress that 
in order to preserve energy conservation within our model, 
there are a number of possibilities regarding the kinetic en- 
ergy within the ejected wind. As stated above, such possibil- 
ities include that the wind remains in the AD vicinity and 
therefore does not require excess kinetic energy for escape 
from the BH potential. Another possibility is that the wind's 
excess kinetic energy is liberated from photons which are not 
along our line of sight and therefore does not affect the ob- 
served SED. Including all such considerations, there are two 
main results of the new calculations. The first is the under- 
standing that disc winds are able to considerably soften the 
optical-UV SED of accretion discs. Second, we have shown 
that the estimation of Mbh cannot rely solely on measur- 
ment of luminosity, as L/L e dd can differ substantially from 
Min /M Bd d- 

The SED of a wind ejecting disc is influenced by the 
fact that removal of hot accreting gas from the inner regions 
of the disc removes energy mainly from the UV part of the 
SED and has less effect on the optical and NIR continua. 
The extent of this energy removal depends on M; n , M ou t 



and on the shape of the wind. M ou t represents the accretion 
rate over most of the disc and therefore adjusts the peak 
wavelength of the SED and contributes radiation mainly at 
low frequencies, with a lesser contribution to the optical- 
UV continuum which depends on Mt n and rj/2. Further- 
more, the contribution to the long wavelength continuum is 
almost frequency independent whereas the contribution to 
the UV is frequency dependent. Mi n governs the shape of 
the UV continuum such that low values of Mi n compared 
to M ou t soften the UV SED significantly. r 1( / 2 governs the 
relative importance of accretion rates at the inner and outer 
radii of the disc. Specifically, a large value of ri/2 causes en- 
ergy to be removed at larger radii, influencing the SED over 
a large range of wavelengths from NIR to UV. A small 
causes an almost constant accretion rate throughout almost 
the entire disc and a sharp fall of accretion at small radii, 
influencing mainly the UV SED. These factors together, de- 
termine the spectral shape of a wind ejecting disc and are 
evident in figures [4] - [5] We find that we can model the ob- 
served a ° « —0.5 for Ms = 1 BHs and a range of 

1450-4200A ° 

M in , M out , a and r l / 2 - Such soft spectra cannot be mod- 
elled by standard ev-discs except for very large BH mass or 
very low accretion rates. In particular, some of our models 
(1.8 a nd 2.7) fit relatively well with the lVanden Berk et alj 
l|200ll ) composite shown in fig. [S] We note that these mod- 
els have extremely high ratios of M ou t/Mi n and that models 
1.5 and 2.5 (which are standard a-disc models with the same 
Min) already have relatively soft optical-UV spectra. 

Since eqn. [8] is not valid for wind ejecting discs, Mbh 
cannot be evaluated solely by measurement of L. As fig. 
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[7] shows, there is a range of values, Mi„/M e dd, for a given 
L/Ledd- We have found that estimating L by using a certain 
XL\ and a bolometric correction may result in inaccuracy 
by a fact or ~ 2. Furthermore, we find that the method sug- 
gested bv lDavis fc Laorl l^OHl ) for estimation of M by mea- 
suring L v at optical wavelengths is not consistent with our 
wind ejecting models and may cause significant discrepen- 
cies in evaluation of Mbh- Together with knowledge of L, 
information on the spectral shape of the AGN can further 
constrain Mbh- In fig. comparing the left and right pan- 
els, one sees that the slope, a, in some spectral ranges is more 
sensitive to Mi„ than in other spectral ranges. For example, 
a at UV wavelengths is far more sensitive to changes in ac- 
cretion rate at the ISCO, than at NIR wavelengths. This is 
true for all values of Mg, a and M ou t- Therefore, for known 
Mx and a values, measurement of a ° or a 

456-912A 912-1450A 

can remove much of the uncertainty in Mbh caused by dif- 
ferent wind profiles and leaves only the uncertainty caused 
by the unknown value of M ou t- 

Finally, this work is not intended to explore the physics 
of real AGN winds or their influence on the disc geometry 
and physical state. As explained above, there are various 
types of winds (radiation driven, thermally driven, MHD) 
and a large range of factors which may influence the out- 
flows. Our work shows that such winds affect the disc lumi- 
nosity, the SED and the BH growth rate, under very general 
assumptions about the actual wind profile. The understand- 
ing of the effects of disc winds on the SED and on our abil- 
ity to measure L allows one to constrain the possible errors 
made in estimation of BH growth rate. 
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